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ABSTRACT 

Evolved mass-losing stars such as Mira enrich the interstellar medium (ISM) significantly 
by their dust-rich molecular wind. When these stars move fast enough relative to the ISM, 
the interaction between the wind and ISM generates the structure known as the astropause (a 
stellar analog of the heliopause), which is a cometary stellar wind cavity bounded by the contact 
discontinuity surface between the wind and ISM. Far-infrared observations of Mira spatially 
resolve the structure of its astropause for the first time, distinguishing the contact surface between 
Mira's wind and the ISM and the termination shock due to Mira's wind colliding with the ISM. 
The physical size of the astropause and the estimated speed of the termination shock suggest 
the age of the astropause to be about 40, 000 yr, confirming a theoretical prediction of the shock 
re-establishment time after Mira has entered the Local Bubble. 
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Introduction 



Mira is the archetype of asymptotic giant branch (AGB) stars, which contribute si gnificantly to th e 
enrichment of the interstellar medium (ISM) by a dust-rich molecular stellar wind (e.g. ISedlmavii 11994 ). 
Hence, tracing mass- loss histories of these stars has a profound impact on knowledge of the composition of the 
ISM and its evolution in host galaxies. Mira is actually a wind-accreting binary system (jReimers fc Cassatella 
19851 ): a mass- losing AGB star is the primary star (Mira A), dominating the system. Mira's wind is, therefore, 
a cool molecular wind from the p rimary, in which dust grains convert radiation pr essure of the star into the 
outward momentum of the wind i Kwok 1975t Steffen. Szczerba. fc Schhbernerlll998) . 



Previous investigations in CO line emission revealed that Mira's wind at ~ 4 km s^^ formed a circumbi- 
nary molecular envelope of 0.3 in radius, whose otherwise STp herical structure was punctuated by bipolar 
outflows at ~ 8 km s~^ (jjosselin et al.ll2000l : iFong et al.l 120061 ). Mira's space velocit y is approximately 150 
km s ~^ nearly due south as cal culated via a spherical trigonometric tra nsformation ( Johnson fc Soderblom 
19871) using its proper m otion (jPerrvman et al.lll997l ). radial velocity (jEvand 119671) and the solar motion 
lOehnen fc Binnevlfligj ). 



Hence, Mira's wind has been interacting with the ISM flow local to Mira for quite some time. This 
i nteraction left a tu rbulent wake of appr oximately 2 degrees l ong, observed in emission in the ultraviolet 
([Martin et al.l 120071 ) and 21 cm H I line (jMatthews et al.l 120081 ). The ultraviolet emission is thought to be 



- 2 - 



due to H2 excited collisionally by hot electrons arising from the shocked ISM flow. The 21 cm hne emission 
of atomic hydrogen is, therefore, expected from dissociation of H2 into H. 

This is the second case of the wind-ISM interaction around an AGB star afte r R Hya, in which a bow 
shock resulting from the interaction was seen in the far-infrared ( Ueta et al. 2006h . While the H2/H I wake 
of Mira has been attributed to the wind-ISM interaction owing to Mira's large space velocity, distinct shock 
layers of the interaction themselves hav e not been confirme d observationally, except for the possible bow 
shock structure seen in the ultraviolet (^M artin et al. 2007 ). Thus, it has been unclear how the II2/II I 
turbulent wake seen at the large spat ial scale is physicall y related to the CO wind seen at the small spatial 
scale, except that br ight H2 streams ([Martin et al.l 120071 ) are correlated spatially with CO bipolar outflows 
|josselin et al.ll200d l. 



In this letter, we present Mira's far-infrared image taken with the Spitzer Space Telescope {Spitzer; 



Werner et al 



2004), which resolv es spatially the structure of Mira's astropause (a stellar analog of the 



heliopause. iNeutsch fc Fahij 119821 ) for the first time. Resolved are (1) Mira's circumbinary dust envelope 
at the core of the astropause, (2) the termination shock where the stellar wind is slowed down by the ISM 
and (3) the astrosheath in which a turbulent wake assumes the characteristic cometary shape beyond the 
termination shock. Thus, the Spitzer 160/im data, in conjunction with previous multi-wavelength results, 
provide the most comprehensive view of the wind-ISM interaction for an AGB star to date. 



Observations and Data Reduction 



Mira happened to be located at the nort hern edge of the XM M-LSS field in the Spitzer Wide-Area 
Infrared Extragalactic Legacy Survey (SWIRE; iLonsdale et aTll2003l ) . Far-infrared mapping of Mira and its 



vicinity was perfo rmed on 2004 August 1 and 2 with the Multiband Imaging Photometer for Spitzer (MIPS; 



Rieke et al.ll2004|) . Due to the relative position offset among the MIPS detectors in the focal plane, Mira was 
scanned over only in the 24 and 160/Ltm bands (scans 08 and 09; AOR KEYs 5844992, 5845248, 11774976 and 
11775232). These scans were done at the medium rate along five scan legs of 3° long each with a cross-scan 
stepping of 148". Quick inspection of the archived data indicated that Mira was saturated in both of the 24 
and IdOfim bands: however, the 160/^m data appeared salvageable. 

Data reduction was done in the following steps with the Ge Reprocessing Tools (GeRT; ver. 20060415 
of S14 processing]^ and the the Mosaicker software (ver. 20070615)0 First, we used GeRT to create our own 
basic calibrated data (BCD). Bright sources can corrupt stimulato r-response calibrati on if they happen to 
be on the detector when the calibrating stimulator is illuminated (iGordon et al. 2005h . Thus, we excluded 



stimulator frames affected by bright Mira to optimize the time-dependent detector responsivity calibration 
for our particular data set. We also adjusted the number of reads in the data ramp to recover valid data 
in otherwise saturated pixels. Then, we applied a high-pass time median filter with a 50-detector-pixel 
width on the custom-made BCDs to remove the "streakings" due to residual slow-response variations of the 
detector, while masking out a circular region around Mira of 30-detector-pixel radius to avoid filtering out 
the extended structure around Mira. Finally, the custom-processed BCDs were mosaicked into a single map 
with the Mosaicker. 



^Available from |http://ssc. spitzer . caltech . edu/mips/gert / 
^Available from |http://ssc. spitzer. caltech. edu/postbcd/| 
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The MIPS 160/im band is known to suffer from a short-wavel ength light leak, wh ich can result in a 
false point source offset from the position of the true I60/im source iGordon et allboosh . Inspection of the 
mosaicked Mira map indicated that the emission core was elongated by the presence of a secondary point 
source offset from Mira (Figure [T]a,) . Since Mira's J magnitude is —0.7 (more than 6 magnitudes above the 
limit over which the leak becomes detectable), the core elongation was most likely due to the leak. Hence, we 
removed the light leak in the Mira map by scaling and subtractingthe image of a calibration K2III star, HD 
163588, whose 160/im emission was known to be mostly the lealo. This yielded a single-peaked, PSF-like 
emission core. This method would inevitably subtract the photospheric component of HD 163588 (estimated 
to be 56 mJy) from Mira. However, it was deemed negligible with respect to the flux of Mira at 160/im (52 

Jy)- 



3. Results and Discussion 

Figure [T] displays versions of the Spitzer/MIPS image of Mira at 160/im in log-scaled false-color. These 
maps were made at the 12'.'0 pix~^ scale (11.'7 x 13.'3 field of view) with the mean sky coverage of 3.5 ± 1.9 
pixel" ^. The effect of saturation was not removed entirely by adjusting reads in the data ramp, and the 
12'.'0 pix~^ scale turned out to be the minimum pixel scale for which all pixels register valid data. The one tr 
sensitivity and the average large-spatial-scale sky emission (the component removed during the median filter) 
were found to be 0.8 and 8.8 ± 1.0 MJy sr"-'^, respectively. The derived sky emission value was consistent 
with the estimated value of 6.5 MJy sr which is obtained with the Spitzer Planning Observations Tool 
(SPOT) software0 Contours are 80, 60, 40, 20, 10, 5, 2.5 and 0.5% of the emission peak, with the lowest 
contour corresponding to 3 fi (= 2.5 MJy sr^^). 

The light-leak-corrected image is presented in Figure [Uj, in which the emission core is evidently single- 
peaked. The leak-corrected emission profile of Mira can be dissected into three distinct emission regions: the 
emission core, plateau and halo. The core is slightly extended with the full-width at half-maximum (FWHM) 
of 1^0 X 0'.9 at the position angle (PA; measured East from North) of 70°. Since nearby point sources have 
the average FWHM of Of 7, the extension of the core appears genuine. The plateau is an elliptical region 
(3f4 X 2'.8 at the PA of 70°) of relatively flat profile encircling the core (at 5% to 30% of the peak emission). 
The halo is a fainter nebulosity surrounding the plateau. If we define the periphery of the halo at 3 cr (the 
lowest contour in Figure [T]), then its extent is 8.'5 x 6.'5 at the PA of 15°. The scan-map direction is the PA 
of 160°, along which some residual "streakings" are still recognizable below 3 a. Hence, the elongation of 
the 160/im emission around Mira is not the artifact of scan map ping. Mira has bee n imaged in the far-IR 



most recently with the AKARI Infrared Astronomy Satellite ((Murakami et al.ll2007[ ). AKARFs Mira maps 



are consistent with the Spitzer map with a single-peaked emission core with extension at the PA of 15° (H. 
Izumiura et al. 2008, private communication). 

The point-spread-function (PSF) in the MIPS 160/im band has the first Airy pattern out to ^ 7' (at 
12" pixel"^) at ~ 4% of the peak surface brightness. A synthesized PSjH was convolved with an elliptical 
Gaussian function to conform the observed shape of the core and subtracted from the data (Figure [T];). We 
immediately see that there is still substantial emission (~ 50 MJy sr^^ 2Qa) in the plateau, suggesting 



^MIPS Data Handbook Ver. 3.3.1 (|http: / /ssc. spitzer. caltech.edu/mips/dh/[ ), p. 76. 
*http: / /ssc. spitzer. caltech.edu/documents/background/ 



^Available at http: //ssc. spitzer. caltech.edu/archanaly/contributed/stinytim/ 
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that it is caused by an ellipt ical rinK of 2.4 x 2'.0 at the PA of 70° . The core is oriented in the sam e way as the 



circumbinary CO envelope ijjosselin et al.ll2000), while the CO en velope is smaller than the core (jFong et al 
[2006) . At the distance of 107 pc to Mira ( Knapp et al. 20031 ). the physical size of the CO envelope is 
3 X 10^^ cm in radius. This agrees well with the expected ra,dius of Mira's molecular envelope resulted from 



mass loss at the rate of - 2 x lO"^ Mq yr'^ |josselin et al.lboool: iRvde fc SchdieilboOll: IPong et allbood) 



including the effect of photodissociation by the interstellar radiation field ( Mamon. Glassgold. fc Hugginj 

Therefore, the core represents Mira's circumbinary dust envelope whose physical size is ~ 4.5 x 10 
cm in radius. 

Since the Mira system is moving at ^ 150 km s~^ th e resulting wind-ISM interaction has given rise 
to the ultraviolet bow shock structure (jMartin et al.ll2007l l. despite Mira's slow wind velocity of 4 km s~^ 



|josselin et al.lboool: iFong et allbood ). Figures [2^ and [2l3 clearly show that the downstream side of the 



ultraviolet bow shock coincides with the southern edge of the far-infrared halo. Moreover, the cometary shape 
of the halo is the most consistent with the bow shock interpretation, in which the far-infrared nebulosity 
is shaped by a turbulent wake flowing almost due north. According to the theory of the (solar)wind-ISM 
interaction, the ISM flow approaches supersonically toward Mira, becomes subsonic past the bow shock 
and streams around the astropause, while Mira's wind expands freely to form a bubble bounded by the 
termination shock, beyo nd which the wind becomes compressed and turbulent and flows downstream in the 
astrosheath (Zank 1999[ ). Assuming that the 160/im emission probes Mira's dusty wind emanating from the 
circumbinary envelope, we interpret that the leak-corrected, core-subtracted image distinguishes spatially 
all components of the astropause (Figure [H:) with the halo tracing the astrosheath and the emission ring 
delineating the termination shock (see, also. Figure [2ji). 

The prominent ultraviolet streams (Fig ures and[^b) are outflows of H2 excited collisionally by hot 
electrons in the bow-shock-excited ISM flow jMartin et al. 2007). These H2 flows, therefore, must have pen- 
etrated the termination shock and been directed downstream while getting e xcited in the astrosheath. In the 
astro sheath, H2 can be dissocia ted by either the interstellar radiation fiel d ( Mamon. Glassgold. fc Hugginsl 



19881 ) or collisional excitation (jMartin et al.l 120071 : iMatthews et al.l 120081 ). The observed H I emission is 



mostly concentrated in the vicinity of the termination shock, with the peak FWHM of l.'S ([Matthews et al 
[2OO8), which is compatible with the size of the molecular radius given the resolution at 21 cm (Figure [2j;). 
Similar to II2 , atomic H is directed toward downstream beyond the termination shock, lending strong support 
for the multi-phase characteristics of the flow in the astrosheath. 

The stagnation point of the wind-ISM interaction is therefore at the apex of the halo (astropause), 
which is 3.6 x 10^^ cm from Mira. Requiring ram pressure balance at the stagnation point, this implies that 
the IS M density local to M ira is nisM = 0.02 cm~^. This value is consistent with the result of a two- wind 
model (jWareing et al.ll2007l ) that accounts for Mira's entry i nto the Local Bu bble, the region of warm 10® 



K), tenuous (~ 0.01 cm'^) iSM in the solar neighborhood (lLallementll200ll ). Given Nm < 1.3 x 10^^ cm 



and 25% abundance of neutral species ( Matthews et al. 2008[) . we can estimate the optical depth at 160/im 
to be 8 X 10~® assuming a power- index scaling law with the index of 1.2 (i.e. silicate dust). Under the 
assumption that thermal dust emission dominates at 160/im, the emission map and the optical depth yield 
the temperature of dust; ^ 500 — 1100 K in the circumbinary dust envelope, ^ 200 K in the termination 
shock and ~ 30 — 200 K in the astrosheath. 

Based on these estimates, the pre-shock temperature of Mira's wind is thought to be ~ 50 K (assuming 
gas-d ust thermalization), implying Mach number of 7.5. With the generalized jump relations (e.g. iPetelskil 
19811 ). the velocity of the termination shock is —2 km s^^ in Mira's frame and the ratio of post- to pre-shock 
densities 3.8. These values yield the post-shock gas temperature of ~ 900 K. Thus, in the shock-excited 
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astrosheat h, dust is still present (the gas temperature is below the condensation temperature of silicates 
« 1200 K: ISalpeteiill977l ). but probably de-thermalized from gas. Spectroscopic follow-ups are necessary to 
investigate the evolution of the dust-gas characteristics upon passage of the termination shock. 

At 2 km s^^, the termination shock should have taken ~ 40, 000 yr to advance to the present position 
(1 X 10^^ cm) from the astropause (3.6 x 10^^ cm). The total flux of the astrosheath (including the termination 
shock) is 37 Jy, and this translates to the total mass of 7 x 10~^ Mq using the dust emissivity of generic 



interstellar dust JPraine fc 1^2007*). With Mira's mass loss rate (-- 2 x lO""^ Mq vr~^: Ijosselin et al.llioooi 
Fong et allbooeHRvde fc Schoier.200l M it should have taken ~ 40, 000 yr for the astrosheath to form. These 
estimates for the age of the termination shock is consistent with the time needed for the termi nation shock to 
re-est ablish itself after entering the Local Bubble (^ 40, 000 yr) as estimate d numerically bv IWareing et al 



(j2007l ). who attributed a kink in Mira's ultraviolet tail ([Martin et al.ll2007^ to Mira's entry into the Local 
Bubble. 

The presence of the bow shock in dicates that the o uter shock is supersonic. In the warm (~ 10^ K), 



tenuous (~ 0.01 cm '^) Local Bubble (|Lallementl 120011 ) ■ however, the shock velocity of ^ 150 km s ^ is 



barely supersonic. This suggests that the pre- and post-shock parameters across the bow shock are nearly 
continuous. This is consistent with the fact that the ultraviolet emission profile does not indicate any density 
enhancement (i.e. material pile-up) at the southern edge of the bow. Hence, ultraviolet emission in front of 
the upstream face of the astropause is due to excitation of the post-shock ISM gas by a shear flow. While 
Mira's wind and the ISM gas flow are in general separated by the astropause, these flows are inherently 
turbulent. Thus, these flows are mixed in the downstream wake and II2 is coUisionally excited, resulting in 
the observed far-ultraviolet emission. Excited II2 can be dissociated, leading to the observed H I emission. 
Dust grains help H2 to replenish in the wake, and this cool, turbulent three-phase gas flow results in emission 
of II2 and H I continuing at least for 1.5° long. 

Hence, there is an emerging picture of Mira's astropause (Figure [2li) . The extent of the far- infrared 
emission outlines the astropause, which is the contact surface between Mira's molecular wind and the ISM 
flow. The bulk of the halo emission is due to dusty molecular wind material in the astrosheath, shock-excited 
by the termination shock delineated by the emission ring. The termination shock is the discontinuous inter- 
face between the pre- and post-shock regions of Mira's wind, inside which Mira's wind expands undisturbed 
coming off the circumbinary envelope hosting bipolar outflows. 

Since the post-shock regions of both Mira's wind and the ISM flow have a significant physical size, the 
shock must be non-radiative (i.e., it has not been cooled down and compressed into an unresolvable size). 



According to iFerland et al.l (jl998l ). in a coUisionally excited solar wind abundance gas heating dominates 
cooling below lO'* K (their Fig. 6). The estimated post-shock gas temperature of Mira's wind is ~ 900 K. 
Thus, the p ost-shock gas can indeed maintain its t emperature, i.e., its physical size. Because previous numer- 
ical models ( Wareing et al ] l2007l : iRaga et al. I l2008h were concerned mainly with the structure of the turbulent 
wake of the wind-ISM interaction, follow-up multi-phase hydrodynamical modeling at low temperatures to 
reproduce the structure of the astropause appears extremely worthwhile. 

The total flux of the astropause i s 52 Jy, which im plies the total mass of 1 x 10^^ Mq using the dust 
emissivity of g eneric interstellar dust (jPraine fc Lill2007i) . Since the total atomic mass in the astropause is 
1.3 x 10^^ Mq ( Matthews et al.llioOS ). the total molecular mass is estimate to be 8.7 x 10^^ M© to be roughly 



"Various authors quote a range of mass-loss rates for Mira, which is a variable star. In the present paper, however, we are 
concerned with the long-term effects of Mira's mass loss. Therefore, we adopt the average rate of mass loss in our discussion. 
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1 to 9 atomic-to-molecular mass ratio. The observed far-ultraviolet luminosity implies the H2 dissociatio n 
rate of ^ 2.5 x 10^^ s"""^, if collisional excitation accounts for all far-ultraviolet emission iMartin et al."2007'). 
This rate, however, would have yielded a factor of four more atomic H than observed (.Matthews et a l. 2008). 
Dust grains can potentially provide H2 formation sites to reduce excess atomic H. Using parameters of gas 
in the astrosheath (nn « 7 and Tgas ~ 900 K) and sihcate grains (0.1/xm radius, 3.7 g cc~^, 1% mass ratio), 
the H2 formation rate is approximately 4 x 10^" s~^, if one assumes that an encounter with two atomic H by 
a dust grain always results in production of H2 . This suggests that H2 far- ultraviolet emi ssion is not entirely 



due to collision of H9 with hot electrons and by some other means as already implied (jMartin et al.l 12007 



Matthews et al.ll2008r ) and/or there exist some mechanisms to bolster production of H2 from atomic H. 



Clearly, Mira's astropause presents a unique laboratory for investigations of a multi-phase hydrody- 
namical flow, turbulence and dust-gas processing, which has a profound impact on the chemical makeup 
of th e ISM. Because recent findings suggest that wind-ISM interaction around AGB stars may be com- 
mon (jUeta et al. 20061 2007), spatially resolved spectroscopy of multi- phase dusty, low- temperature gas flow 



will be extremely productive in the coming er a of Herschcl, SOFIA and ALMA (e.g. Pilbratt et al. 2003t 
Becklin. Tielens. fc Cal"li]l2007t loiofssonlboosl . respectively). 



Ueta thanks T. Lc Bcrtrc and L. D. Matthews for making their H I data available and R. E. Stencel 
for careful reading of the manuscript. This work is based on archival data obtained with the Spitzer Space 
Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology under a 
contract with NASA. Support for this work was provided partially by an award issued to the University of 
Denver by JPL/Caltech. 
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Fig. 1. — Spitzer/MIPS far- infrared image of Mira's astropause at 160/xm. (a) Calibration-optimized, 
saturation-corrected image in the 12"pixel~^ scale () centered at the position of Mira (north is up, east to 
the left) based on SWIRE/XMM-LSS scans, which yielded 1 a sensitivity of 0.8 MJy sr'^ The contours 
represent 80, 60, 40, 20, 10, 5, 2.5 and 0.5% of the peak surface brightness, which is indicated at the bottom 
right corner in MJy sr~^. The lowest contour corresponds to emission at 3 cr level (= 2.5 MJy sr~^). The 
extent of the detected (> 3a) emission halo is 8.'5 long and 6'.5 wide roughly at the PA of 15°. (b) Same 
as (a), but for the short-wavelength light-leak removed image. The slightly-elongated, singly-peaked core 
(> 40% of the peak) is surrounded by the similarly elongated plateau (> 5% of the peak) at the PA of 70°. 
(c) Same as (b), but for the core-fitted-PSF subtracted image. The lowest contour is 2.5% of the peak. The 
plateau is not entirely due to PSF and there is a ring of emission (at ~ 40% of the peak; 2.'4 x 2.'0 at the 
PA of 70°). 




RA offset (arcmin) 



Fig. 2. — H2/H I imaging data overlaid with the 160/ini contours of Mira's ast ropause, and a sch ematic of 
its structure emerging from the present study, (a) Near-ultraviolet color image ([Martin et al.ll2007n overlaid 
with the leak-corrected 160/im contours, showing the spatial relationship between emission in these two 
bands. The contours represent 80, 60, 40, 20, 10, and 0.5 % (3 a) of the peak surface brightness, (b) Same 
as (a), but far- ultraviolet color image (jMartin et al.ll2007l ) ov erlaid with the leak-co rrected, core-subtracted 
160/Ltm contours, (c) Same as (a), but 21 cm H I color image ( Matthews et al. 20081 ) overlaid with the leak- 
corrected, core-subtracted 160/J,m contours, (d) Schematic of the structure of Mira's astropause revealed by 
the 160/im image and its spatial relationship with respect to H2/H I distribution. The circumbinary dust 
envelope, the termination shock and the astrosheath are spatially distinguished in the far-infrared for the 
first time. 



